Natural versus forced convection in laminar starting plumes 
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A starting plume or jet has a well-defined, evolving head that is driven through the surrounding 
quiescent fluid by a localized flux of either buoyancy or momentum, or both. We studied the 
scaling and morphology of starting plumes produced by a constant flux of buoyant fluid from a 
small, submerged outlet. The plumes were laminar and spanned a wide range of plume Richardson 
numbers Ri. Ri is the dimensionless ratio of the buoyancy forces to inertial effects, and thus our 
measurements crossed over the transition between buoyancy- driven plumes and momentum- driven 
jets. We found that the ascent velocity of the plume, nondimensionalized by Ri, exhibits a power 
law relationship with Re, the Reynolds number of the injected fluid in the outlet pipe. We also 
found that as the threshold between buoyancy-driven and momentum-driven flow was crossed, two 
distinct types of plume head morphologies exist: confined heads, produced in the Ri > 1 regime, and 
dispersed heads, which are found in the Ri < 1 regime. Head dispersal is caused by a breakdown 
of overturning motion in the head, and a local Kelvin-Helmholtz instability on the exterior of the 
plume. 



INTRODUCTION 

Jets and plumes are flow structures of considerable in- 
terest [1] |2] due to their widespread occurrence in in- 
dustrial and natural systems, from fuel injection [3 , to 
mantle convection [H |5l El |7] • Both jets and plumes occur 
when flow discharges from an isolated, submerged source. 
The distinction between jets and plumes is that a pure jet 
is driven only by momentum flux at the source, while a 
pure plume is driven only by buoyancy. If a discharge has 
a combination of both momentum and buoyancy, there 
is no sharp distinction between a jet and a plume, and 
flows span a continuum of possibilities between the two. 
Such a flow structure may equally well be referred to as 
a "buoyant jet" or a "forced plume" ; we adopt the latter 
terminology in this paper. We experimentally studied 
the morphology of laminar forced plumes and establish 
the general scaling of their ascent velocity over the full 
range spanning pure jets and plumes. We also identify a 
sharp transition in the shape of the plume head as the 
type of forcing is varied. This transition can be traced 
to the onset of a divergent flow structure in the under- 
belly of the head. Surprisingly, this transition does not 
interrupt the scaling of the ascent velocity. 

The degree to which buoyancy affects flow emerging 
from a source is usually quantified by the Richardson 
number Ri, defined in Eq. [T] which is the ratio of buoy- 
ancy forces to inertial effects. Thus, Ri can also be 
viewed as a measure of the extent to which convection 
is "forced", in the Ri <C 1 limit where inertial effects 
dominate, as opposed to "natural" or "free" convection 
which occurs in the purely buoyant, Ri 1 limit. In our 
experiments, slightly buoyant (< 1% density difference) 
fluid was injected into an ambient, quiescent fluid of con- 
stant density. A wide range of Richardson numbers were 
accessed by varying the fluids and the flux of buoyant 
fluid at the source. We measured the ascent velocity and 



morphology of forced plumes across this range. 

Not only do jets and plumes occur in a wide variety 
of phenomena at various scales, they are also studied in 
a vast array of geometric configurations and flow scenar- 
ios [El [9]. We focus here on unconfined flows, not di- 
rectly interacting with a boundary, and on flows taking 
place in a quiescent ambient medium of uniform density. 
By contrast, other important cases are jets and plumes 
in a crossflow or in a stratified fluid [9 . For a given 
set of experimental conditions, the most important fac- 
tor determining jet or plume behaviour and evolution is 
whether the flow is laminar or turbulent. In addition to 
being either laminar or turbulent, jets and plumes can 
be a mixture of both types of flow, with near-field lam- 
inar flow that gives way to a far-field turbulent flow [8 . 
We focused on laminar, forced, compositionally buoyant 
plumes, for which precise control and characterization of 
morphology is possible. 

Forced plumes can be positively buoyant, in the sense 
that buoyancy forces act in the same direction as the in- 
jected momentum, as we consider in this paper, or they 
can be negatively or neutrally buoyant, as in the pure 
jet limit. More complex forced plumes can have several 
sources of buoyancy, such as in double-diffusive salt fin- 
gers [10^ or chemically reacting plumes [TTJ [I2l [13] . In 
this paper we only consider buoyancy produced by com- 
positional differences. Finally, the state of jet or plume 
evolution involves one more important distinction: be- 
tween a steady, well-developed flow and a starting plume, 
the transient case for which the injected fluid is penetrat- 
ing the ambient medium and growing in length. Over the 
course of its evolution, a laminar starting plume develops 
its most visually striking feature: a well-defined, evolving 
head. 

Given the wide variety of types of jets and plumes, 
and the wide array of scientific contexts in which they are 
studied, it is no surprise that there is a lack of uniformity 



in the terminology used for jet and plume anatomy. As 
shown in Fig. [TJ we adopt the term head to describe the 
structure on top of the plume, and conduit to describe 
the fluid emanating from the source that connects and 
provides volume flux to the head. Some authors refer to 
a starting jet or plume head as a cap, and the conduit is 
sometimes called the stem^ or corridor in the literature. 
The conduit starts at the outlet of a small pipe, through 
which buoyant fluid is being pushed. Finally, we use the 
term lobe to describe the bottom, under-turning part of 
the head. 

The most prominent feature of a laminar starting 
plume is the vortex ring that often forms in the plume 
head. The generation and evolution of vortex rings has 
been a subject of longstanding interest in fluid dynamics 
due to their natural beauty, their utility in engineering 
applications and their rich history as simple solutions of 
the fluid equations [M]. Vortex rings can be produced 
experimentally by using a cylindrical piston to inject a fi- 
nite volume of neutrally buoyant fluid into quiescent sur- 
roundings [15]. Continuously supplied, neutrally buoy- 
ant jets have also been used to study the velocity fields 
of evolving vortex rings [16j . In addition to neutral buoy- 
ancy scenarios, vortex ring formation and pinch-off has 
also been investigated in the context of buoyant starting 
plumes, where results suggest some universality between 
free vortex rings produced by pinch-off and vortex rings 
generated by a piston [171 HH] • Buoyant vortex rings may 
also be created by the pinch-off of accelerating plumes 
driven by chemical reactions [11 . 

The majority of work on laminar jets and plumes has 
focused on the steady variety. Nonetheless, starting jets 
have attracted interest in studies of combustion [19] , and 
starting plumes in the context of geophysical applica- 
tions [11 [5j m I3- Previous work on laminar starting 
plume heads has focused on ascent velocity [20l [211 IIS] , 
temperature and concentration measurements [23], ve- 
locity field measurements [24 , and scaling laws for the 
head [21]. Since plumes in these studies were driven by 
a thermal buoyancy flux produced by a localized heater, 
a relevant dimensionless quantity is the Prandtl number 
Pr= V / ti^ where v is the kinematic viscosity and n is the 
thermal diffusivity. The dependence of plume flow on Pr 
has been investigated [22 . For plumes driven by compo- 
sitional differences, the molecular diffusivity D replaces 

and the relevant ratio is known as the Schmidt num- 
ber, Sc= v/D. Typically, D <^ so that the diffusion 
of compositional differences is very slow on the timescale 
of the evolution of the plume. This is the case in our 
experiments, for which 2 x 10^ < Sc < 2 x 10^. 

Diverse plume experiments have yielded a number of 
plume head morphologies, such as cam^?/ structures [Sll^, 
umbrellas [7 , heads that encapsulate a vortex ring, like 
the one shown in Fig. [l} and heads that do not contain 
vortex motion in their lobe [21]. Some thermal plume 
heads become non-axisymmetric under sufficiently high 
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FIG. 1: The anatomy of a laminar plume. The entire plume 
is shown on the left, while a magnified image of the head 
is shown on the right, h is the height of the plume from 
the outlet to the top of the head, Wh is the width of the 
head, and ih is the head length. The plume shown is from 
set D5 and was created using an injected flow rate of Q = 
1.33 X 10~^mL/s. In this image, h — 19.3 cm, Wh — 1-9 cm, 
and ih — 1.6 cm. 



forcing [25 . More recently, a numerical model of ther- 
mally driven 2D line plumes found four different dynam- 
ical regions of starting plume morphology [26j. 

While a variety of head structures are known, the na- 
ture of transitions in the type of head that forms as ex- 
perimental parameters are varied has not been systemat- 
ically explored. Moreover, the conditions under which a 
stable laminar head changes its morphology have received 
little attention. Similarly, the scaling of the ascent ve- 
locity of forced buoyant plumes has not been previously 
studied. We address these issues in this paper. 

In the next section, we describe the experimental ap- 
paratus and protocols, and in the succeeding section, we 
present the results for the scaling and head morphologies. 
This is followed by a brief conclusion. 
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EXPERIMENT 



The apparatus used for our forced plume experiments 
was a plexiglass tank with a vertical glass capillary tube 
built into the centre of the tank floor. The inner diam- 
eter of the capillary tube, which served as the outlet for 
fluid injected into the tank, was 3.0 mm, and the in- 
ner dimensions of the square tank were 13.4 cm between 
walls that were 50.2 cm high. A second, smaller square 
tank was also used which had 9.8 cm between its inner 
walls and a height of 33.7 cm. A syringe pump was con- 
nected to the bottom end of the outlet pipe, allowing fluid 
to be injected at a controlled, steady rate. The syringe 
was filled with glycerol-water solutions slightly less dense 
than the ambient glycerol-water solution into which they 
were injected. Ambient solutions used in this experiment 
ranged from a 20% to an 80% volumetric ratio of glyc- 
erol/ water. The densities pa and pi of the ambient and 
injected solutions were measured using an Anton-Paar 
densitometer. Density measurements were made at the 
temperature that had been recorded in the room dur- 
ing the respective experiments. The room temperature 
was measured to ±0.5° C, implying a possible system- 
atic error in the absolute densities pa and pi of less than 
±0.0003 g cm~^. However, the difference in densities, 
= Pa — Pii is rather insensitive to the temperature, 
because of the very small difference in the thermal ex- 
pansion coefficients of the two fluids. Since the ambient 
and injected fluids were nearly identical, the dimension- 
less ratio of their kinematic viscosities, Vijva^ was typi- 
cally close to unity. The relevant dimensionless density 
difference /S.p/ pa^ along with other properties of the ex- 
perimental fluids, are given in Table |lj 

Two methods of visualization were employed to ob- 
serve the plumes. In the smaller of the two tanks, shad- 
owgraphy was used. This technique requires a constant 
light source to be directed horizontally at the tank. By 
virtue of the difference in the refractive index of the in- 
jected and ambient fluids, the less dense (injected) fluid 
projects a dark image onto translucent white tracing pa- 
per attached to the tank at the opposite side to the light 
source. A CCD camera was used to capture the shad- 
owgram images of the ascending plume. The sets of ex- 
periments for which shadowgraphy was used to observe 
plume behaviour are specified in Table [l] by set names 
beginning with "S". In the larger tank, blue dye was 
added to the injected fluid for visualization purposes. As 
fluid was being injected and a growing plume was formed, 
a digital CCD camera was used to capture images of its 
evolution. These are given in Table |l] by set names begin- 
ning with "D" . For both the dye and the shadowgraphy 
experiments, the injected flow rate Q was varied in the 
range of 3.3 x 10~^mL/s to 6.67 x 10~^mL/s. Altogether, 
34 forced plumes were analyzed. 
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FIG. 2: The vertical position of the plume head as a function 
of time for the D4 experimental runs at five injection flow 
rates. The solid lines show the linear best fit for plume heights 
above 4 cm. The slopes of these lines are the constant ascent 
velocity vu of each plume. 



RESULTS 
Ascent velocity and scaling 

After a short transient, plume heads ascend at a con- 
stant velocity. Linear fits of h values extracted from time 
lapse images were used to determine head velocities, 
as shown in Fig. [2] These fits excluded data in the im- 
mediate vicinity of the outlet, which was selected to be 
h < A cm. Typically, in the very early stages of plume 
head formation, the ascent of the head is slower than in 
the linear regime as plume accelerates towards the con- 
stant velocity that it eventually achieves. 

As is apparent from Fig. [2j for the D4 set of experi- 
mental runs, the ascent velocity of the plume head, 
increases with the injection flow rate Q. This trend is 
consistent for all data sets. To establish the scaling of Vh 
for all experimental runs, we applied dimensional analy- 
sis. We assumed that the flow in the system was depen- 
dent on five physical variables: <i, ^i, Q, and g' . Here, 
d is the diameter of the outlet tube, Vi is the kinematic 
viscosity of the injected fluid, Q is the injected volume 
flux, and g' = gAp/pa is the reduced gravity, where g 
is the acceleration due to gravity. Since there are five 
variables and three fundamental dimensions, a straight- 
forward application of the Buckingham 11 theorem [27 
implies that two dimensionless groups describe the sys- 
tem. Choosing some convenient numerical factors, the 
dimensionless groups can be written 
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Here, Ri is the plume Richardson number, the ratio of 
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Set 


pi 


pa 


Ap/pa X 10^ 


Vi X 10"^ 


l^a X 10-^ 






(g/cm^) 


(g/cm^) 




(mVs) 


(mVs) 




SI 


1.0590 


1.0599 


0.0849 


1.80 


1.86 


0.968 


S2 


1.1151 


1.1166 


0.134 


3.95 


4.17 


0.947 


Dl 


1.1745 


1.1775 


0.255 


16.5 


17.8 


0.927 


D2 


1.1752 


1.1797 


0.381 


16.5 


19.1 


0.864 


D3 


1.1897 


1.1920 


0.193 


24.3 


26.3 


0.924 


D4 


1.2045 


1.2160 


0.946 


40.7 


66.2 


0.615 


D5 


1.2135 


1.2161 


0.214 


59.9 


66.2 


0.905 



TABLE I: Fluid properties of the various injected (subscript i) and ambient (subscript a) glycerol-water mixtures. Densities 
were measured with a densitometer at the same temperature at which a set of experimental runs were performed. The viscosities 
were interpolated from data given in Ref. [28j . 



the buoyancy forces driving the plume to the inertial 
terms in the Navier-Stokes equations. Re is the stan- 
dard Reynolds number for the flow in the outlet pipe 
based on its diameter d and average flow velocity Vavg- 
We have Q = Avavg in a pipe with the cross sectional 
area A = (7i/A)dP. 



10" r 




FIG. 3: The dependence of the plume Richardson number on 
the Reynolds number of the injected fluid flow in the outlet 
pipe. The solid curve is a power law fit to the data, which 
gives Ri = 4.3 Re~°"^^. The dashed red line indicates Ri=l. 
The symbol shading indicates the plume head morphology; 
shaded symbols indicate confined heads, while open symbols 
denote dispersed heads. 

Plotting Ri vs. Re on logarithmic axes, as shown in 
Fig. |3j we find that all the experimental data collapse on 
a single curve. The dependence is well described by a 
power-law that spans almost three decades in Re more 
than two decades in Ri. Re covers a range from well 
inside the Stokes regime, Re ~ 0.1 to Re ~ 100, but 
remains small enough that the flow in the outlet pipe 
is never turbulent. The range of Ri straddles Ri ~ 1, so 
that we observe both flow regimes that are dominated by 
injected momentum (Ri <C 1, for large Re), and plumes 



in which buoyancy forces dominate (Ri 1, for small 
Re). The power law is of the form 



Ri = a Re^ 



(2) 



From a least squares fit of all of the experimental data 
on Fig. [3) it was determined that a = 4.3 ± 0.2 and 
k = —0.96 ± 0.05. Thus, we arrive empirically at the 
remarkably simple result that Ri is approximately pro- 
portional to Re~^, or equivalently that Ri Re = const. 
From Eqn. [l] this implies 

=4.3 ±0.2 , (3) 



Ri Re 



i.e. that all forced compositional plumes are described 
by a single dimensionless group that is independent of (i, 
the diameter of the outlet pipe. 
Isolating Vh in Eqn. 3, we find 



vy, = (0.54 ±0.01) 



1/2 



(4) 



The above expression for the head ascent velocity may 
be regarded as a generalization for forced compositional 
plumes of the classic scaling proposed by Batchelor [29] 
for the velocity Vc of purely thermal plume conduits given 
by 



\vpCp 



1/2 



(5) 



Here, the thermal buoyancy flux gaP/ pCp replaces the 
densimetric compositional buoyancy flux g'Q in Eqn. [4j 
where a is the thermal expansion coefficient, P is the 
power input by the heater, v is the kinematic viscosity 
of the isoviscous solution, p is the density of the ambient 
solution, and Cp is the specific heat of the fluid. For the 
centreline velocity of a steady thermal plume, Vcu tfie 
proportionality constant in Eqn. [5] is a known function of 
the Prandtl number |30j. 
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a) 

FIG. 4: The nearly self-similar evolution of a typical confined plume head. Each image is 5s apart. The plume is from set D5 
and the injection rate was Q — 2.00 x 10~^ mL/s. From the first to the last image in the sequence, h increases from 20.5 cm 
to 35.2 cm, and the size of the head increases from 1.9cmx 2.3cm {^h x Wh) to 2.4cm x 3.1cm. 



where e is a solution of e^loge"^ = Pr~^. This result 
could presumably be generalized to the forced composi- 
tional case for which the Schmidt number Sc would re- 
place Pr. 

Experiments [21] [22] on thermal starting plumes in flu- 
ids with various Prandtl numbers found that the head rise 
velocity Vh was related to the conduit centerline velocity 
given in Eqn. |6]by 



(0.57±0.02)vc^. 



(7) 



The prefactor in Eqn. [7] is less than one for the obvious 
reason that the growing plume head must be supplied by 
a higher speed conduit. 

The physical reason for the scaling of the ascent veloc- 
ity in each case is clearly the fact that the morphology of 
the plume eventually becomes independent of conditions 
near the isolated source of buoyancy flux. In our case, 
this implies that should become independent of the 
diameter of the outlet pipe d. Since d is the only length 
scale in the problem, only a single dimensionless group is 
required in the overall scaling of v^^ as in Eqn. [4j 



Morphology of the plume head 

The simplicity we found in the previous section for 
the scaling of the ascent velocity does not extend to the 
scaling of the plume head morphology. While the as- 
cent velocity shows the same power-law scaling across a 
wide range of Re, the head morphology cannot be de- 
scribed by a single characteristic scaling over the same 
range. We found instead that there were two distinct 
types of plume head that exist on either side of Ri ^ 1. 
For lack of an established taxonomy, we classify these 
head morphologies simply as confined for Ri > 1 and 
dispersed for Ri < 1. The distinguishing feature is the 
presence or absence of a stable vortex ring in the lobe of 
the plume head as it grows during its ascent. Confined 
heads, which are observed for larger Ri, exhibit such sta- 
ble overturning structures, while for unconfined heads, 
observed for smaller Ri, the vortex ring exists for only a 
short time before it collapses and disperses. The domain 



12 






1 r 






• 


3.3x10"^ mL/s 


10 




□ 


6.7x10"^ 






♦ 


1.33x10"'' 


8 




A 


2.67x10"'' 






T 


4.00x10"'' 



0) 



6 - 



4 - 



2 - 




J- 



4 6 



10 



FIG. 5: Head width as a function of head length for the D4 
set of experiments. All of the plume heads in this set were 
confined. The head scales wu and lu have been nondimen- 
sionalized as Reynolds numbers Yleuw and Re^^, respectively. 



of each type nicely meet at Ri ~ 1, as indicated in Fig.jSj 
which forms a reasonably sharp boundary between the 
two morphologies. In addition to the stability and size 
of the vortex ring generated in the lobe, the two types 
of head have other distinguishing characteristics that are 
described below. 

A confined head is the classic mushroom-shaped lami- 
nar plume head. This type of head is shown in Fig.[l] and 
in a sequence of time-lapsed images in Fig. [4] from a D5 
experiment. We use the term "confined" to describe this 
type of head because the fluid that comprises it remains 
within a compact structure. This structure is preserved 
for the duration of the life of the starting plume, which 
ends when the head collides with the top of the tank. The 
fluid in the head circulates around an axisymmetric vor- 
tex ring that remains localized near the top of the plume. 
As shown in the sequence in Fig. [4j in order to accommo- 
date the influx of new fluid delivered to the head through 
the conduit over time, a confined plume head grows and 
its vortex ring increases in size. We studied the growth 
of the head by measuring the width of the head and 




FIG. 6: A sequence of images of a dispersed plume head during the evolution of a D2 starting plume with Q — 2.67 x 10~^ mL/s. 
Each image is 2s apart. From the first to the last image in the sequence, h increases from 10.4 cm to 31.0 cm. The size of the 
plume head in a) is = 1.6 cm and Wh — 2.1 cm. 



the length of the head as defined in Fig. [T] at var- 
ious times. These two lengths were nondimensionalized 
by their associated Reynolds numbers, 



Vh^h 



and Re 



hi 



(8) 



Figure [5] shows how the Rehw scales with Re^^, for all 
of the confined D4 plumes. We find that Rehw is simply 
proportional to Re^^ independent of Q. In all cases, 
R^hw ^ C ReMi with C = 1.24 ± 0.04, and hence the 
aspect ratio of the head is constant, with Wh/^h = 1.24 ± 
0.04. This behavior is typical for plume heads in the 
confined regime. 

This simple scaling is not found in dispersed plume 
heads. The evolution of a dispersed head is shown in 
Fig. [6j In contrast to confined heads, dispersed heads do 
not remain compact and do not contain a stable vortex 
ring structure. Instead, the height of a dispersed head 
elongates faster than its width, as shown in Fig. [7| for 
Dl plumes. The head dimensions Wh and ^h are only 
clearly defined in the early stages of growth, after which 
instabilities in the lobe take over, as shown in Fig. [6j 

A vortex ring structure forms early in the evolution of 
a dispersed plume head, but subsequently becomes un- 
stable by an interesting mechanism. The vortex in the 
lobe fails to draw in all of the fluid being delivered to 
it from the outer layer of the head. Instead, the flow 
diverges and some of the fluid is directed upward and 
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FIG. 7: Head width as a function of head length for the Dl 
set of experiments, nondimensionalized as Reynolds numbers. 
With the exception of the plume with the lowest value of Q, 
all these plumes are dispersed. For all dispersed plumes, the 
head length grows faster than the width as the lobe becomes 
unstable. 



eventually escapes the vortex ring. At the moment when 
this vortex entrainment breaks down a hammer shaped 
structure, unique to dispersed heads, is observed. This 
structure is barely discernible in Fig. [6k,b, and is shown 
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in detail in Fig. |8] At the site of the hammer struc- 
ture, the plume head protrudes outwards, forming a bell 
shape. Above the protrusion, an instability, almost cer- 
tainly of the Kelvin-Helmholtz type, develops in the outer 
fluid layers, and the vortex ring becomes cut off from the 
influx of new fluid. Thereafter, the remnants of the vor- 
tex ring near the lobe dissipate, and an elongating fluid 
skirt that trails below the advancing head is formed. The 
evolving skirt exhibits further Kelvin-Helmholtz type in- 
stabilities, while the overall flow remains axisymmetric, 
as seen in Fig. [6] A contributing factor to the devel- 
opment these instabilities are small amplitude bulges of 
fluid that develop in the plume conduit. The bulges de- 
velop well above the outlet, and move upward at a higher 
velocity than v^^ causing fluctuations in the volume flux 
feeding the head. Bulges are visible in the portion of the 
conduit shown in Fig. [6]f-j. 




FIG. 8: An image of the axisymmetric hammer-shaped struc- 
ture (indicated by the arrows) that results from the onset of a 
divergent flow structure beneath the underbelly of the head. 
This structure is not observed in confined heads. The head is 
from a D3 plume with injection rate Q = 6.67 x 10~^ mL/s. 
The dimensions of the plume are h = 21.9 cm, Wh = 3.7 cm, 
and Ih — 3.3 cm. 

The scenario for dispersed head evolution we have de- 
scribed above is typical of heads that form near Ri 1. 
For Ri <C 1, the head forms Kelvin-Helmholtz instabili- 
ties relatively quickly, before the hammer and bell shapes 
have time to develop. It should be noted that a Kelvin- 
Helmholtz instability will develop when the local shear 
reaches Ri < 1/4 [T, where the shear velocity and thick- 
ness of the shear layer are the important parameters. It 
is known from simulations that below this critical 
value of Ri, Kelvin-Helmholtz instabilities appear in a 
starting plume head. The global Richardson number we 
have defined for the whole plume, Eqn.[T] is not identical 
to the local Richardson number that governs the Kelvin- 



Helmholtz instabilities of the plume head. 



CONCLUSION 

In summary, we have experimentally explored the scal- 
ing and morphology of forced compositional plumes, 
which could also be called buoyant jets, in the lami- 
nar regime. We focused on starting plumes, for which 
a well-defined, ascending head exists. From dimensional 
analysis, we found that the system is described by the 
Richardson number of the plume and the Reynolds num- 
ber of the injected buoyant fiuid in the outlet pipe. We 
experimentally determined that the Richardson number, 
which scales with the inverse square of the head ascent 
velocity, had a simple power law relationship with the 
Reynolds number, which scales the volumetric fiux of the 
injected fluid. This scaling is exactly such that the di- 
ameter of the outlet pipe d drops out. This reflects the 
physical fact that the ascending compositional plume's 
morphology and speed eventually become independent of 
the details of the localized source of buoyancy and mo- 
mentum that produced it. This result generalizes some 
previous observations and theory for thermal plumes to 
the case of forced compositional plumes. Our results are 
speciflc to the case of large Schmidt number and nearly 
isoviscous plumes. 

The morphology of the advancing plume head exhibits 
two clear forms, depending on the Richardson number. 
For large Richardson number, we observe conflned plume 
heads which contain a stable vortex ring and retain 
their mushroom shape throughout their evolution. These 
plume heads show a simple self-similar scaling. For small 
Richardson number, the plume heads become unconflned 
when the flow feeding the vortex ring fails to close and a 
thin, trailing skirt is formed. This skirt undergoes sub- 
sequent local Kelvin-Helmholtz instabilities. 

It would be interesting to extend this study to the case 
of forced plumes with both thermal and compositional 
effects, to non-isoviscous cases and to forced plumes as- 
cending in a density stratifled medium. 
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